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ON CERTAIN DETERMINANT FORMS AND THEIR APPLICATIONS.* 
(First Paper. ) 


By Pror. W. H. Ecuoxs, Charlottesville, Va. 


Determinants no longer occupy in analysis the position of mere symbols 
of tabulated results. In the higher analysis they enter largely into the opera- 
tions themselves. Not only are they powerful levers in these far-reaching 
methods of investigation, but are as well the instruments for discovering new 
relations which may exist between quantities and functions when bound to- 
gether by their side lines. While there may be little or nothing that is new 
in the results here obtained, the operations serve somewhat to illustrate in an 
elementary manner the sweeping comprehensiveness of determinant forms. 

Consider the alternant determinant function, 


Je, 1,2, #, 


2 
Tan, 1, Any + 


(1) 


2 
Jay, 1, yy Toy a," 


Let the function 7% and its first x derived functions be finite and contin- 
uous for all values of the argument « from a, to a,,. 


Let 
< & < i 


< 
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be quantities such that 
< 


< 


I. GENERALIZATION OF ROLLE’s THEOREM. 


The modern form of Rolle’s theorem is this: “If a continuous function 
fr vanish when # — a, and also when a == 4, then its derived function 7’2, if 
also continuous, must vanish for some value of 7 between a and b.” The proot 
of it goes hand in hand with the conception of a derived function. 

Lagrange’s form of Rolle’s theorem is this: “If #2 be a continuous fune- 
tion for all values of between 7 and — then 


where w is some value of « between @ and 6.” 

It is proposed to generalize this form of the theorem, and to express 7" 
in terms of Fins Gg +++» 

The determinant function /% vanishes whenever a takes any one of the 
n 4 1 values «,,...,¢,. Its first derived function, therefore, vanishes when- 
ever z takes any one of the ~ values ,..., 4, —;. So on, until finally its nth 
derived function vanishes for some value w of 7, such that a, << u < a,. 

The nth derived function of F’ is, as may easily be seen, 


Fe = + (—1)"t'n! 


where J/ is the minor of 77, and J/, that of 2" in the determinant /’r. 
Hence 


n—1 | 
=n! 
1, Any ’ a, ’ ar 
n 
1, Qs » MW 


which is the desired relation. 
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If we represent the minors of fw,,..., fa, in the determinant J/, by 


A,, ..., Ay, respectively, we have 
' 
fm =F (A, fa, + A, fa +... +4, (2) 


where M and the A’s contain only the a’s, and are therefore independent of 
the form of the function. 


I.—If fa, = fa, =... — fa, = fa, 

then 
Ke — M (fe + fa), 
and 
Fre = Mf'e; 

whence the derivatives of % and fx vanish together. 

Since 

Fz=0, for @,,...,5 Gq 

we have 


F%, —0, FG, =0,..., Fb_, =9, 
In like manner, 


= 0, 
The same results hold in the particular case when 
which is the extended form of Rolle’s theorem due to M. Ossian Bonnet. 


IT.—If all of the different values of fu, except fu,, fa,, fit, ..., 
are equal to each other and to fu, we have 


+ A, fa, + A, fa, + A, fa, 
= A, (fa, — fa) + A, (fa, — fa) + Ar( fu — fa), 
since (A, + A, + ... + A,) is equal to the determinant formed by replacing 


each of the fa’s by unity, and therefore vanishes. 
If all are equal except one, fu,, then, writing out the binomial difference- 
products of the alternants,* 
A, n!’ 


= (— 1)" (a, — a) — (4, — ay) * 


* Muir’s Determinants, p. 162, et seq. 


Sima 


108 ECHOLS. ON CERTAIN DETERMINANT FORMS AND THEIR APPLICATIONS. 


Corresponding results hold when the fv values vanish as a particular case 
Thus, putting x for a,, and 


Ja, Ja, =O, ..., =% 


we have 


Sa = — a,) (@ — a)... —a,)° 


a form also due to M. Ossian Bonnet. 


When x = 1, we have the following geometrical illustration of Lagrange’s 
form of Rolle’s theorem ; 
| fo, 
Fx 


Fay, a, 1 


is the double area of the triangle whose base is the chord (fa, a@,), (fa, a) of 
the curve fx, the vertex of the triangle being a running point on the curve. 

This area vanishes when the vertex coincides with either extremity of the 
chord. So that, if the function fr is finite, continuous, and single-valued 
between a — a, and 7 —a,, we have the first derivative of Fr vanishing for 
some value ~, of x, between a, and a,. Hence, 


Ja, — fay = (a, — %) f’u. 


The chord being regarded as of constant length, the change of Fr is 
therefore proportional to the altitude of the triangle. As the length of the 
chord converges to zero, the altitude of the triangle for z = « converges to the 
mid-ordinate of the segment of the osculating circle on that chord. 


IL. 


If all of the » | 1 values of fa are known except fu,, we have by (2) 


Writing out the alternants into binomial difference-products and cancelling 
common factors in numerator and denominator of the coefficients, we have 


— ay) (a, — (4 — Ay) (a,—@y) (4, 
(a,—d)) (a,—q)) (a,—4n_1) 
(4,—) 


+ — (a — (a, — a) 


n! 


— 
| 
n! 
| 
j — 
| 
| 
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If we leave off the last term, which renders the equation exact, we have 
Lagrange’s interpolation formula, the most important one in analysis, the 
principle of which is to substitute for the unknown function a rational integral 
function of the (x — 1)th degree which has in common with it the 2 values 
Jay, -++, fan. The residual term in f'”w corresponds to the remainder after 
n terms in Taylor’s formula, and shows that if the unknown function be a 
rational integral function of no higher degree than x — 1, Lagrange’s formula 
is exact, since f"~ — 0. If the unknown function be of the wth degree and 
the coefficient of x" be unity, then the formula above gives exact results, since 
f*e= ni. 

The above result was written down merely to identify (1) with Lagrange’s 
formula. The interpolation formula is written out at once from the original 
determinant ; thus, letting J/,, M,,..., If, be the minors of 1, z,..., ”" in 
Fr, we have 


(3) 


In this, / vanishing for 7 = “,,..., @,, makes fx and the rational inte- 
gral function ¢v, of the nth degree, have » ~ 1 common values. Neglecting 
Fx/M, we interpolate any function fx for a given x between a, and a,. Com- 
puting the coefficients in gx once for all, any number of such interpolations 
may be rapidly made. 

Fx/M being a function which vanishes for the » + 1 values 7 =a, ... 
a,, We may write it 


dy) (@ — a)... — ay) 


The form of ¢/a* depending in general on that of 7. 

If fr be a rational integral function of the nth degree, so also must 7x 
be. But Fz vanishes x -- 1 times; it is therefore zero for all values of .r, and 
the interpolation is exact. If fr be of the (m -+ 1)th degree, so is Fr. Hence 
yx cannot contain z, and is constant. 

If in (2) we put the factors in the coefficients of fa,, etc. equal, so that 


a, — a, = 4,--a,=—... = 
we have 
An f™u fa San 


*It is shown below that 7 = (— 1)"+1 f 


n+l), 
(m+1)!’ where <0 < 


o , 
niga M 
| 
on | 
| 
J M 
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If A/n = 1, then 


= fa, — nfo, fay + (— Wn, 


a formula analogous to the fundamental interpolation formula of finite differ- 
ences. 


In (3) the coefticient V,/M of «” is the same as that of a", if in M, we 
change a, into ,, and multiply by (— 1)"~”. 


ITT. MecuHanicaL QUADRATURE. 


The area of the curve fx is, by (3), 


a, a 


fe de = ede + 
Jo 


ay ay 


Approximately 


dz = dz , 


the error being 
f Ix 
M 


The curve cuts the axis +1 times at the points == a, ..., dy, 
and the area is made up of positive and negative portions which tend to annul 
each other. It seems probable that with a proper distribution of the n — 2 
intermediate points a, ..., ,, and frequently of but one or more of them, 
the error area may be made to vanish. 

Gauss ( Werke, Vol. III, p. 203) shows how to make this distribution inde- 
pendently of the form of fz, when 7x is a rational integral function of a degree 
not higher than 2a (see Jacobi’s proof; Boole’s Finite Differences, p. 52). 
Gauss’s quadrature formula is the integration of Lagrange’s interpolation ~ 
formula, in which the error area is 


fc 1)" (a, — a)... (a, — ay) da,. 


ay 


i 
| 
| 
| 
| 
| 
| 
| 
| | 
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As a simple illustration of the error vanishing when fr is an irrational 
function, the approximate formula gives exactly the quadrant of a circle 
(radius = 1) for a, = 0, 4, = 1, and a, having some value between 37 and 
#73, the error for these values being -} 0.0006 and — 0.00165, respectively. 

Of course, if the ordinates be equidistant, the general formula reduces to 
the ordinary well known forms. 

The most important researches in the theory of interpolation have had 
reference to this formula. The investigations of Minding, Christoffel, Mehler, 
and Jacobi (all published in Crelle’s Journal; see Boole, Finite Diff., p. 56) 
have had in view the general problem of determining the position of the ordi- 
nates we should choose, so that the approximation may be as close as possible. 
There still remains much to be done in this direction. The great importance 
of the computation of definite integrals is sufficient cause for close research in 
the attempt to reduce the error area for the general definite integral. 

We may examine 
a, 


in the following manner :— 
The original determinant suggests that we may put 


1 


2 n+l 


Fe JR, 


2 n+l 
1, dy, Gy 


When « has some fixed definite 


where /2 is some unknown function of . 
value 7, we have 


Fx, 
R= 
which is independent of x 
The function 
Fx — Sh, 


vanishes when 2 has any one of the 2 values 4,,...,@,, Its 4. 1)th 
derivative 


Fz — 


0 dx n+l J, 


— +1)! MR,, 
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will therefore vanish for some value of 7, say v, such that v lies between a,, and 
a,, the greatest and least of these values.* Hence 


(n+ 1)!" 
Since ~, may be any value of « between a, and a,, we may drop the suffix, and 
put 


* It is easy to show that Prof. Sylvester's general alternant 


S12 S 2%, 


may be treated in the same manner, obtaining the relation 


| 
} 
| 


in which the last constituent of each column except the rth and the sth is zero, and 


n+l 
u 

gutly = anti 
u 


where 
n+l n+l n+l n+l 


u 


in which the rth column is wanting, and 


| Sete, seta, ..., ...+ 


in which the sth column is wanting. 
The quantity w has some value between the greatest and least of the values z,, ..., 2%, and 
is to be substituted for z after the (m + 1)th differentiation. 


| 
| 
i 
| 
| 
| 0, 
| 
| 
i ‘ 
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Therefore 


ay a, 
Fr 
n+l 

ay ay 
the same result as before.* 

Any effort to reduce this, when the form of the function is known, must 
be made through an inquiry into the value and form of v,+ which is a function 


of F(x), a, and 
FENERALIZATION OF THE EXPANsion OF fr. 


We begin by deducing Taylor’s formula in determinant notation and then 
proceed to the generalization of the expansion of the function 7x. 

1. Zaylor’s Formula.—It is required to find a relation, if such exists, 
between the two values of a function, 7 and f(« + h), the first » integral 
powers of /, and the first » derivatives of fr. 

Write down in a line, 7( + /) followed by the successive powers of 4 up 
to the nth inclusive. Underneath this line write a second line formed by put- 
ting / = 0 in the first line. Underneath this second line write a third line 
formed by differentiating the first line with respect to /, and in the result put- 
ting A= 0. Repeat this operation until the first line has been differeutiated 
n times. Draw the side lines, thus forming a determinant function J). This 
determinant 7 is a function of the quantities whose inter-relation is desired. 


+h), 1, hy 


fe, a@.... 0 

0,110, 0, ..., 90, 0 

fre, 0,0,2140,..., 0, 0 
D=|*. 

0,0,0, 0, 0 

0,0, 0, 0, .... (m,—I)!, 0 

fx, 0,0,0, 0, ..., 0, n! 


We observe that ) and each of its derivatives down to the nth, inclusive, 
vanishes if A vanishes. 


*See Synopsis der Hoeheren Mathematik von Johann G. Hagen, Erster Band, pp. 158, 207, 
or Baltzer, Theorie der Determinanten, p. 87, where the alternant /'z may be found, and notes on 
Gauss’s quadrature. 

+ See a paper on the remainder after x terms in Taylor’s formula, by A. W. Whitcom, Ameri- 
ean Journal of Math., Vol. III, No. 4, which indirectly bears on this subject. 
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Factoring out 1!,2!,..., !, from the 2nd, 3rd, ... , (2 + 1)th columns, 
respectively, we have* 


I(x h), 1, A, 314", aye” 
Se, 1, 0, O, ..., 0 

0, 0, 0, Cony 1 


Observing that in the determinant the general term is 


Ere, 


we have 


n | D 
This series is convergent if / is finite, since by increasing x sufficiently 
we may make 7)/n!! as small as we choose. 


The form of the determinant suggests putting 


Aer 


in +1)! 


R, 


where /2 is some unknown funetion of « and A. 
When / has some definite fixed value /,, we have 


or 


which is independent of A. 
Then the function 


+ h), 1, A, nih”, ame 


& & 0 
Se, 0 
Sve 0 
R., | 1 


*T use the symbol n!! to represent the double factorial product 12!3!...n! 


| 4 
if 
| 
ine 
| | 
| 
im 
| 
| = 
| 
Ant! 
i 0 
| | 11 
| 
aq 
| | 
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and each of its derivatives to the nth, inclusive, vanish when 2=—0. But 
the function vanishes when / == /,. Its first derivative vanishes therefore for 
some value /, between 0 and /,. Its second derivative vanishes also for A,, 
such that 0 < A, </,; and so on, until finally its (2 —- 1)th derivative vanishes 
for some value of A, say such that 0 — 4h —h, < hy. 

The (x — 1)th derivative of the above determinant is 


kh) 
hence 

R, (a +- th). 
We have, therefore, 


D, 
0 in + 1)!" 


Oh). 


Dropping the suffix, because /, is any finite value of A such that fe and 
its first (x + 1) derivatives are finite and continuous between the limits 2 and 
x -- h, we have, finally, 

Je, 0 
0h),0,0,...,0, 1 


The original determinant was filled in with zeros merely in order to sim- 
plify its expansion. It might, however, have been written in the more general 


form 
1, (@ +A), (@ + Ap, ..., (@ 


Je, lL, a, 
Sj*2, 0, 0, 0, 


and the result in corresponding shape. 

2. Generalization of the Expansion of fre.—While, in substance, the preced- 
ing method is but Maclaurin’s theorem, yet it may be regarded as an extension 
of it, inasmuch as it affords quite an independent method of expanding any 
function in powers of a variable. For while, after Maclaurin’s method, we 
have successively after differentiation put the variable equal to zero, we might 
have put it equal to any constant, in order to make the successive derivatives 
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independent of the variable. Or, indeed, we may substitute for the variable, 
after differentiation, any series of constants,* as in the following :— 

Let gx be any function of ~, finite and continuous, together with its first 
n derivatives, for the values a,,..., a, of the variable. 

Write down ¢ followed by the » powers of 7. Underneath write these 
values when « =a, Underneath write the derivative of the first line, and 
substitute in it a, for #, and so on until the nth derivative is written with «, 
for 2. Draw the side lines and call the determinant / 

Factor out the factorials in the diagonal as before, giving the result ; thus, 


gv, 1, a, x" 


1 2 1 n 


Let 
1 1 +1 
1.2 1 n 1 n+l 


way & Uy 
(—1)"*' al! JR. 
Where 7? is some unknown function of ur. 
When « has some fixed finite value «,, we have 
Ry =(— 1 
which is independent of 
The function 
F—(—1)""' al! 
or 
ge, ..., 


1 n+l 


ge, 0, 1, «.. 


n! 


or in general any functions, which are independent of 2. 


i 
| 
i @ 
| 
| 
i 
i 
| 
| 
i. Throughout this paper I mean by constants, in the sense as employed above, any quantities, 
| 
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vanishes when « = «,, also when w — «,; therefore its first derivative vanishes 
for some value of 7, say 7,, between a, and 2,. Since its first derivative van- 
ishes for + =, and also for # = a,, its second derivative must vanish for 
some value of 2, say ”,, between «, and a, And so on, until finally its 
(n + 1)th derivative vanishes for some value of .7, say u, between .”,, and a,. 


Hence 
Fy —(— 0. 


In this equation we may drop the suffix, since #, is but avy fixed finite 
value which « may have, within the limits of continuity. We have, then, in 
general, 

F—(—1)"' all OO, 
or 
gz, 1, 2,..-, 
ga, 1, 


ie ei 0. (5) 


We bring ¢’a, to the first place in the first row by 7 | 1 interchanges ; so 
we have for the corresponding term in the expansion of the determinant 


and when transposed it becomes 
(— r) A, ¢"a,. 
Whence the series 


gx = ga, — A,¢'a, + A,g’a,—...+ (—1)" A, ¢"a, + (—1)"" (6) 


in which the coefficients A,,..., A,4,, are the minors of ga, ..., ¢"4,, in the 
above determinant, and are independent of the form of the function. 

If all of the values «,, @,,..., 4, become equal to any one of them, say 
a, or if the values @,, a, ..., 4, become equal to x, the series becomes Tay- 
lor’s expansion; and if a, = a, == ...@, = 0, it becomes Maclaurin’s. — If 
dy 0, and a, = a, =... =a, x, we have John Bernoulli’s series. 


— 
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In order to avoid the tedium of expanding the determinant coefficients 
A,,..., A,, we observe the following :— 


1 2 1 

1 

Qa, My a, 


2 1 1 1 
1 r-l 1 r—1 | 
1 r-2 1 r—2 | 
q 0, 0, ’ 0, 0, 1 
1 
1 
\ 0, 1, 5) de. 
a, | 
0, 0, 0, 1, a, 
But 
1 2 1 r—l | 
1, My > 
1 
1 r—2 
» 


0,0, 0 ae 


Therefore, we observe that any A may be derived from the preceding one 
by an operation analogous to that of integration. Integrate the variable ac- 
cording to the regular rules for integration, and at the same time increase all 
the suffixes of the arbitrary constants by one. Thus, using the symbol I to 
denote the operation,* we have 


A,=— A,_, de. 


* The symbol I and that of the inverse operator gy should in addition include the change of 
algebraic sign. 


| 
| 
| 
| 
0,0, 0, at | 
| 
| 
| 
| 
i 
| 
| 
| 
an 
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The functions ¢’a,, independent of 7, may also be included in the opera- 
tion by increasing both the suffix of the argument and the index of the function 
by unity. Thus 


A,¢’a, = — | A,_1 ¢" 'a,_, dz. 


(—1y | de. 


We write the series, finally, 


= 
t= 


A, dx + g"* u A, de, (7) 
in which the symbol "Id means that the above detined operation is to be 
repeated 1 times. 

A,. is a homogeneous function of the 7 + 1 quantities ,, %,..., @,—,. In 
general A, may be made to vanish by giving either ¢, ,, or as, ... , @ any 
one of the one, or two,..., (7 ~~ 2) values which satisfy A, = 0. For example, 


A, = 0, when a, == a). 


To give the series arithmetical meaning, let « have some fixed definite 
value A; let a, ..., a, be intermediate values in the ascending order 
of magnitude between the fixed limits a, and a, between which limits ¢” and 
its first » derivatives are finite and continuous for all values of 2. 

In A, put first 


the inferior limit. Then A, becomes 


( 1)" (a 


Now put 


the superior limit. Then A, becomes 


By taking r sufficiently large as we may make these two values differ as 
little from zero as we choose. 

In general the value of any co-efficient A, will depend upon the distribu- 
tion of the intermediate values «,,..., a@,,, and may be made as small as we 
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please, or be made to vanish, by letting any one of the a’s have a particular 
value. We have, for computing the A’s, the formula 


1 2 A r 4 


For example, from this we obtain 


A = ay), 


A, = — (# — a,*) + 3, (2? — a”) — a, a, — ay) + (@ — a), 


and so on, 
To investigate the relation which must hold between the a’s in order that 
the series may be convergent, after some fixed term, we have by the above 


formula 

é | r—2 . r—3 r-4 | 

A, 1 3 3.4 d, J 

3 4 
r—2 3 r—3 3 4 r—4 
— 


In the ratio of the second term of the member on the right the first terms 
of the numerator and denominator are identical, while each succeeding term 
of the numerator is less than the corresponding term of the denominator, since 
generally a, >a, ,, the ratio in question is therefore less than unity. 

Letting 4,/A,,-- », we have, therefore, 


also, & priori, 


Therefore, if @, lies between the limits 0 and — 1, or 0 and — 1, the 
series is convergent. 
If we arrange the series according to the powers of x, we have 
peti 


it 
| 
| | 
| 
| 4,4 > > G4 — $ 4,23 
| 
| 
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wherein the determinants C,, ..., (,_;, are independent of 7. We have 


1 2 1 n+l 


4, @y 


as, 0, 1, 1 ay 1 


1 1 2 1 3 
— gay aga, + ay 91% a; ’ + ete 
& 1, a, 
a 


” n+1 9 1 


with the relation 


As before, 
eats O>e>1; 
and 


Again, we have 


C=—¢a,4+ 1,),a7 | 1, | 4 
0, «4, 10, a, 
| 0, 1, a, 
x a x 
Cr = — Cde; | 
201 1 
0 2! 0 
at 


= 


Hence the series* may be written 


<n dr = (—)’ Ce. 


* Since the only condition in general imposed on the quantities @,,..., @,, is that they shall 
be functions independent of 2, a number of series may be deduced from the general one through 
suppositions regarding these arbitraries. While I have only noticed the simple particular cases of 
Taylor’s, Maclaurin’s, and Bernouilli’s formule, I have prepared for publication a paper dealing 
with this series alone in which other forms are noticed, the consideration of which would be beyond 
the design of the present paper. 
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In the determinant (4) let gx be a rational integral function of the 
(n + 1)th degree. And let 


G4 = 0, = 0; 


so that a, is a root of ga = 0, and a,..., a, are roots of its first n derived 
functions respectively. 
Then 
= (n + 1)!; 
and we have the identity 


1 n+l 
1, gem” 

1 n+l 


nt 


If, therefore, we know a root of any equation and a root of each of its 
derivatives, we may at once write down the function from the above. 

Also, since we can always factor out at sight 7 — a, from the right hand 
member of the identity, we can write down immediately the first depressed 
equation of « —= 0, which has for its roots the other x roots of gv = 0, with- 
out actually performing division. 

If —- 0, then 


1 

(0) (n + 1)! 
¢(0) ( 1) r 1)! 0, 1, @2 


is a homogeneous relation between the roots, and the member on the right is 
the product of all the roots of ga = 0, with their signs changed. 

The most general power series of a function in terms of its variable, which 
can be conceived is 


fe = aat + + cat +... + 228, 


where 4, 3,7, ... are positive integers, and a, 4, c, ... are functions independent 
of x. 

If the series be convergent, it may extend to an infinite number of terms ; 
otherwise the number of terms must be finite, and instead of a series we have 
a formula of relation. 


| 
ig 
| | 
| 
| 
it 
j 
| 
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To investigate the above relation, or its existence when there are x — 1 
power terms; write down the rows as before and factor out coeflicients and 
factorials. Thus, we have 


pl’ x!’ et 


(a—n)!’ (8—n)!’ — an)! 


Put 
78 as astn 
ws 
(a—n)!’ (3—xn)!? (6 +y4—n)! 
IR. 


Where is some unknown function of 
When z = ~,, we have 
D, = 4, , 


a relation of quantities independent of ~. 
We have, if 
Dd, 

ie 

4, 
0, for and «=> 
=0, and«=a, «@, between and ; 
=0, for and x, between and «, ; 


| 


= 0, for and r=o,, between 7,_, and w,_, ; 
0, for «= uw between «, and 


124 ECHOLS. ON CERTAIN DETERMINANT FORMS AND THEIR APPLICATIONS. 


Letting );'*' denote the result of differentiating x -- 1 times with respect 
to w the determinant and in the result replacing z by u, we have, with like 
meaning for 


or 
D, = Jy 
Drop the suffix because “, is any value of and put 
Dr 
+1 
(—1)"*' Fu. 
We have finally 
a xs astn 
1 
Iu, 0, 0, 0, i 1 


This is the general expression sought, and determines the forms which the 
coetticients a, 4, c, ... must have. 
More generally still, being given any + 3 functions of 7, such as 


form a determinant Fas before with n + 2 of these functions omitting the last 
¢g. Then form a new determinant exactly as before in which however we 
replace the 7th function 4 by the function ¢, which was previously omitted, 
and call this latter determinant J. 


Put 
F= Jk, 
being an unknown function of 
If z= %, 


is independent of ~. 


| | Di) = Dy, 
| | 
| | 
| | 
| 
| 
} 
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The (x + 1)th derivative of the function 


hy 


vanishes for some value of x, say u, which lies in magnitude between the 
greatest and least of the values a, >, c,..., 2. So that 


R 
n+1 
= 
Put 

— rt 


J are 
Then, dropping the suffix, we have 
F= 4, 


Transposing the 7th column to the first in 7, and to the last in J we have 
the relation* 
Oa Yay Ris Ca 


, , 
bs ay | 


n no gn wn 


0, 0, ...,0, 1 


if the functions «a, f, 7, ... and their derivatives are finite and continuous. 

The great comprehensiveness of this formula is such that it would lead us 
altogether beyond the limits of this papert to do more than consider one form 
of a particular instance of the simplest possible case; namely, that of three 
functions @, a, and 3. 


* While I have written out the first row in the variable as functions of 2, and the second row 
as functions in which 2 is replaced by a, a quantity which is independent of 7, and then proceeded 
with the derivatives, we might have inserted previous to differentiation a number of function rows 
in which 2 is replaced by quantities w, y, z etc, independent of 7, thus making the determinant 
even more markedly an alternant-differentiate. This determinant I venture to name a *‘ composite,” 
jnasmuch as it is composed partially of the form of an alternant and partially of a form which is 
that of a Wronskian. With this distinction, however, in the case of the latter; that the quantities 
under the derived functional signs may in general be entirely independent of each other and are 
independent of the corresponding quantity, z, under the primative functional sign. 

+I have prepared a paper upon the application of the composite to the expansion of an arbi- 
trary function in terms of certain transcendental functions, forms of sin 7 and cos x, which 

presents many interesting points. 
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Here we have, changing the notation to more familiar forms, the relation 
ge sy gy egy (i) 
fugy—pufy 


where u« lies between and y. 
In particular, if the function ¢/ is a constant, say c, then we have 


gu fy —fu gy _ gufy — fugy (ii) 
fu — fy tu 
We know that 


Se — fy — yf 
where v lies between « and vy ; therefore 


gu fy — fegy f% fue 


The member on the left is the form of Professor Cayley’s resultant. 
If ¢ is also a constant, say, /, then (iii) reduces to the familiar form of 
Lagrange’s formula, 


Se — fy = fro. 


University oF Viroerta, Dec. 15, 1891. 
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TISSERAND ON THE PRESENT CONDITION OF THE THEORY OF 
THE MOON.* 


While preparing his “ Traité de Mécanique Céleste ” Tisserand has given 
us many valuable notes on subjects in celestial mechanics. One of the most 
interesting of these is this “ Note sur Etat Actuel de la Théorie de la Lune,” 
in which he has followed pretty closely Professor Newcomb’s “ Researches on 
the Motion of the Moon,” but has presented, however, some important additions. 

Although when Hansen’s Tables of the Moon were published it was 
thought that the complicated problem of the Moon’s motion was definitively 
solved, and although these tables seemed to satisfy very exactly the observa- 
tions from 1750 to 1850, the period of accurate observations which Hansen had 
at his disposal to examine, yet soon after 1860 it became evident that the 
tables were in error, and that the errors were increasing. 

The question how Hansen’s Tables satisfied observations before 1750 was 
taken up by Newcomb in his Researches on the Motion of the Moon. All the 
available material was collected, and carefully and ably discussed by him. He 
finds from the eclipses of the moon given in the Almagest of Ptolemy the fol- 
lowing corrections to the mean longitudes of Hansen’s Tables : 


Epoch. Corr. to Tables. 
— 687 
— 381 —27 +5 
— 189 —20 +3 
+ 134 —16+4 
The eclipses observed by the Arabs give 
Epoch. Corr. to Tables. 
850 — + 
927 —16+1.7 
986 —45+1.3 
For the epochs after 1620 Newcomb finds the following corrections : 
1625 + 50” + 13” 1775 
1650 +39 + 5 1800 0+1 
1675 32 + 1 1825 0 +1 
1700 +21 + 1 1850 0 +1 
1725 + 72+ 1 1875 —8 +1 
1750 0+ 1 


*Note sur l’Etat Actuel de la Théorie de la Lune. Par M. F. ‘Tisserand. Bulletin Astrono- 
mique, November, 1891. 


7 
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In endeavoring to account for these quite inadmissible errors Tisserand 
considers the calculation of the solar perturbations, the calculation of the 
inequalities of long period, and the numerical determination of the constants. 
The coefficients of the inequalities arising from the action of the sun were 
determined by Hansen and Delaunay by quite different methods, and were 
shown by Newcomb to be in practical agreement. 

For the two inequalities of long period arising from the action of Venus 
the values adopted by Hansen are 


V, = + 157.34 sin (— g — 167 + 189” + 33° 36’), period of 273 years, 
V,= 4 21 47 sin (89° — 189 + 4° 44), period of 239 years ; 


while the values found theoretically by Delaunay are 


V, = + 16.34 sin (— g — 16 7 + 18 + 35° 16.5), 
V,= + 0”.27 sin (8 g’ — 18 g’ — 41° 48’), 


where g, 7, 7’ are the mean anomalies of the Moon, the Earth, and Venus 
respectively. 

Hansen’s V, was found by theory, while V, was determined by him em- 
pirically so as to satisfy observations from 1750 to 1850. Also, Hansen has 
adopted 12.17 for the secular acceleration s, though the theoretical value de- 
termined by Adams and Delaunay is 6’.18, the former value being supposed 
better to represent the ancient eclipses. 

When Hansen’s V, is added to the tabular corrections they assume the 
following values, failing to satisfy the observations from 1750 to 1850: 


1625 + 33” 1775 21” 
1650 + 18 1800) + 15 
1675 4-15 1825 + 2 
1700, + 16 1850 — 11 
1725 + 16 1875 — 28 


1750) 4+ 19 


Since the presence of the term J’, furnished Hansen with an erroneous 
value of the mean motion, a correction dn is needed; also, corrections ds and 
de should be given to the secular acceleration and the longitude of the epoch. 
But as an examination of the table just presented shows that there is no sys- 
tem of values of dn, ds, and de that will make the corrections vanish, we have 
to assume that the present theory is incapable of representing the observations 
from 1625 to 1875, and all that can be done is to introduce an empirical term, 
such as 


R =A sin at + Becosat, 
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and then seek to determine the quantities dn, ds, de, A, B, and a, so as to sat- 
isfy the observations. 

Newcomb has used for « a value such that the period of /2 is 273 years, 
which is the period of V,, and the table last given furnishes eleven equations 
of condition for determining de, dn, A, and & in terms of ds. Tisserand solves 
the equations four times, using for the period of /? 391, 273, 209, and 170 
years ; and gives the residuals when ds is left undetermined, and when it is put 
equal to — 6”.0; that is, when « has its theoretical value. The residuals cor- 
responding to the period 170 years of /? are larger than the modern observa- 
tions would admit. The corrections to Hansen’s Tables for 1889.0 are now 
calculated with the first three values of the period of /?, and are found to be 


Period of R. Corr. Corr. 
391 years — 20".1 + ds — 31”.1 
273 years — 15.1 + 0.38 —17 4 
209 years — 78+ 0.69 —11.9 
the values in the second column being those for which ds —- — 6”.0. The cor- 


rection — 17”.4 corresponding to the period 273 years agrees exactly with the 
observed correction. Then with this period, using the corresponding values 
of de, dn, A, and B, are calculated corrections ( and C” to Hansen’s Tables 
for the years 1620-1888, both for s=12".17 and for the theoretical value 
s=6".18. Under V are given the corrections to the tables found from obser- 
vations. 


Cc 0 o-C 

1620 50” + 53’ | 46” + 7’ 
30 48 48 46 
40 + 45 i 43 + 44 —¥ 
50 + 42 + 39 + 42 .3 
60 + 38 + 36 + 39 om & 
70 + 34 -- 33 + 35 — 3 
80 + 30 -- 30 30 0 
1690 + 24 26 + 25 +1 
1700 +. 20 21 + 20 +1 
10 + 15 15 + 15 0 
20 +11 9 +11 — 3 
30 + 7 + § + 6 —j 
40 + + + +1 
50 + jf 0 0 0 
60 0 0 — |} +1 
70 — |} 0 — 3 + 32 
80 — § 0 — 3 + 2 
90 — 0 — +1 
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0’ 

1800 0” 0” 0” 0” 

10 0 0 + 1 —j 

20 + 1 0 + 2 — 3 

30 + j 0 + 2 — 3 

40 + ] 0 + 2 — 3 

50 0 0 + 1] —]1 
50 — 0.7 + 0.0 
52 — 0.5 + 13 + 0.2 + 1.1 
54 — 1.0 + 1.4 — 0.3 + 1.7 
56 - 1.4 + 12 — 09 + 2.1 
58 — 19 + 1.9 — 1.4 + 3.3 
60 2.3 + 2.3 — 1.9 + 4,2 
62 — 2.9 + 2.2 — 3.6 + 4.8 
64 — 3.5 ~ O1 — 3.3 + 3.4 
66 — 42 — 323 — 4.1 +18 
68 4.9 4.0 — §.0 + 1.0 
70 — §.7 — §.4 — 59 + 0.5 
72 — 6.5 — 7.5 — 6.9 — 0.6 
74 7.4 — 91 — 8.0 
76 8.3 9.6 — 9.1 — 05 
78 9.3 — 9.0 — 10.2 +13 
80 — — 10.3 — 11.4 + 11 
82 — 113 — 12.6 — 12.7 + 0.1 
84 — 383 — 14.8 — 14.0 — 0.8 
86 — 134 — 15.4 — 15.4 0.0 
1888 — 14.6 — 16.9 — 16.8 — 0.1 


The computed corrections to the tables agree well with the observed, but there 


would seem to be another inequality having a coefficient of 2” or 3”. 


Having recourse to the eclipses of Ptolemy and those of the Arabs for the 
determination of ds, the most probable value is found to be — 5”.87, or s = 6".3, 
very nearly the theoretical value. 

The most important addition made in this paper by Tisserand to the theory 
of the Moon is to point out that with the addition of an empirical term we can 
represent its motion by the theoretical value of the secular acceleration as well 
as by any. Thus, we can find as good places by using pure theory as by hav- 
ing regard to the friction of the tides, and supposing a diminution in the veloc- 
ity of rotation of the Earth and an apparent acceleration of the movement of 
the Moon. Such a change in the Earth’s velocity of rotation would be so small 
that it could hardly be measured, and though sustained by theory we cannot 
tell whether the change really takes place or whether it is counteracted in 
some way that we do not know. 

The chronological eclipses, which would require s to have the value 12”, 
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are rejected by Newcomb and Tisserand, as the accounts of them are vague and 
untrustworthy. 
The inequality of long period, 


R= A sin at + B cos at, 


which is needed to determine the mean motion with precision, and is intro- 
duced as an empirical term, remains to be accounted for theoretically. The 
combinations that can be made of the mean longitudes of the Moon, the Earth, 
and a planet, and those of their nodes and perihelions are very numerous, and 
there may be, besides J, and the small inequality arising from the action of 
Jupiter discovered by Neison, other inequalities having sensible coefficients. 
Also, the value of J, may not be correctly determined. But Tisserand thinks 
that in the end theory will triumph, and the law of gravitation will completely 


explain the Moon’s motion. 
Asapuo Hatt, Jr. 


NOTE. 


The principal formula demonstrated in the article “On certain space and 
surface integrals,” pp. 61-63 of this volume, had already been given by J. 
Somoff. In the more general case of oblique co-ordinates it appears in his 
memoir “ Moyen dexprimer directement en coordonnées curvilignes quelconques, 
orthogonales ou obliques, les parametres différentiels du premier et du second 
ordres et la courbure Mune surface,” p. 14, Memoirs St. Petersburg Academy, 7 
series, Vol. VIIT. (1865), and is repeated in his Theoretical Mechanics with ap- 
plications to several special cases. It is found on p. 13 of the second volume of 
Professor Ziwet’s translation of the work into German, 7heoretische Mechanik, 
Leipzig, Teubner, 1879. The information contained in this note was kindly 
supplied to the author by Professor Ziwet. 

Tuomas S. Fiske. 


JaNnvuARY 16, 1892. 


ON A NINE-POINT CONIC. 
By Dr. Maxime Bocuer, Cambridge, Mass. 


It does not seem to have been noticed that a few well-known facts, when 
properly stated, yield the following direct generalization of the famous nine- 
point circle theorem :— 

Given a triangle ABC and a point Pin its plane, a conie can be drawn 
through the following nine points : 

(1) The middle points of the sides of the triangle ; 

(2) The middle points of the lines joining P to the vertices of the triangle ; 

(3) The points where these last named lines cut the sides of the triangle. «* 

The conic possessing these properties is simply the locus of the centre of 
the conics passing through the four points A, 2, C, P (ef. Salmon’s Conic Sec- 
tions, p. 153, Ex. 3, and p. 302, Ex. 15). 

Moreover, if we notice that the middle points of the lines A B, AC, PB, 
PC form the vertices of a parallelogram inscribed in the above-mentioned nine- 
point conic, it follows, at once, that the lines BC’ and 7A, being parallel respect- 
ively to two sides of this parallelogram, are conjugate chords of the conie. 
Hence, 

Any side of the triangle and the line joining P to the opposite verter form 
a pair of conjugate chords. 

If each of these pairs of conjugate chords consists of two lines perpendic- 
ular te each other, the conic must become a circle. Therefore, 

Tf the point P lies at the intersection of the perpendiculars dropped from 
the vertices of the triangle A BC upon the opposite sides, the nine-point conic will 
become the ordinary nine-point circle. 

The nine-point conic will be an ellipse when /? lies either within the trian- 
gle A BC or in one of the three infinite portions of the plane which can be 
reached from the interior of this triangle by crossing two of its bounding lines. 
When /? lies in any of the three remaining portions of the plane the nine-point 
conic will be an hyperbola, When / lies ov one of the sides (or extended 
sides) of the triangle A BC, we shall have not a true parabola, as we might at 
first sight expect, but a pair of parallel straight lines. When, however, /? is at 
infinity, we shall have a true parabola (the line at infinity also separating those 
portions of the plane corresponding to ellipses from those correspouding to 
hyperbole). Finally, the case when the nine-point conic is an equilateral 
hyperbola is of some interest, as then the point / must lie on the cireumfer- 
ence of the circle circumscribed about the triangle A BC. 


January, 1802. 
ANNALS Maru. VI, 5. 
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A CYLINDER, diameter 24, intersects a sphere, diameter 27, the surface 
of the cylinder passing through the centre of the sphere. 


Required the part 
of the volume of the sphere contained by the cylinder. 
Artemas Martin. 
SOLUTION. 


Taking the origin at the centre of the sphere its rectangular equation is 


(1) 
and that of the cylindric hole is 


2’ + = Qe. 
Also, 


= ff fae dy dz. 


(3) 
Let « = 7 cos ¢, and y =r sin ¢; whence (1) and (2) become 
r = 2b cos ¢, (5) 
and 
= (f dr dg dz. (6) 
The limits of — (a —/*)and 1 (a@ — 7”); of 7, 2b6cos¢ and 0; 
of ¢, - z and 0. 


1 1 


2” 
4 
0 0 


(7) 


*The above solution is for the case a > 2b. The case a < 2h remains to be treated. The 
ease a = 2b is old.—Ep. 


= cos*¢)! dy. || 


134 SOLUTIONS. 
Let ¢ = : x — #; then dy = — dé, cos ¢ = sin 4, and (7) becomes 
1 
2 
0 
1 
9” 
— : fi (a? — 46° dé 
0 
1 
2 
0 
1 
4 
= 3 za® + 3 a fi (1 — & sin*A) dé 
0 


1, 
2 


16 ab? ( dd. 
0 


Putting e = 2b/u we have 


V = 3 (a 
omit. 
: 
When = 9% 3} [Artemas Martin). 
321 


A ciRcLE meets a hypocycloid of class 3 at six finite points. Show that 
the tangents to the hypocycloid at these six points touch a conic. 


[Frank Morley.) 
SOLUTION. 


The equations of the hypocycloid in circular coordinates are 
a/e=2/t—f ) 1 
y/e=2t —1/#y’ (1) 
t being a complex quantity of modulus 1. (See a paper on the Epicycloid, 
American Journal, Vol. XIII, No. 2). 


Any circle is 
ry sy +7 


| 
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Substituting from (1), we have a sextic to determine ¢, and we observe that the 
product of the roots is 1. 
The tangent to (1) at ¢ is 
uz +vyt+1=d, 


where 
u=—it/e(l+f), v=—F/e(1+P). 


Let the line equation of a conic be 
(A, B,C, G, v, 1° = 0. 
Substituting for ~, v in terms of ¢, we have again 


I] denoting a product. This condition ensures that the six tangents touch a 
conic ; and we saw that it holds if the points ¢ are concyclic. 
[Frank Morley. | 
322 
THE are of a limacon is shown in works on the Caleulus to be equivalent 
to the are of a certain ellipse. Show that the double point on the limacon 
corresponds with Fagnani’s point on the ellipse. [W. B. Richards. | 


SOLUTION. 


It is shown in works on the Caleulus that 


i, 


2 
S=2 ‘(a bY cos*’¢ 4+- (a — bY dg, 
0 


is the quadrant of the ellipse on semi-axes 2(a@ + 4) and 2(a — 6), and also half 
of the whole length of the limacon 7 = a cos 4 - 5. 

It is well known that Fagnani’s point divides the first are into parts whose 
difference is 44, while the half difference between the two loops of the limacon 
is also 40. [W. Echols. | 


cos $0 


Norre.— George R. Dean also points out that at Fagnani’s point cos ¢ = \ = = 
at the node of the limacon.—Ep. 
323 
For solution see Analyst, Vol. I, No. 1, pp. 8-9. I proposed the problem 
in the Schoolday Visitor Magazine for May, 1872, nearly two years before Mr. 
Siverly used it in the Analyst. I had forgotten these facts when I sent it for 
publication in the ANNALS. [Artemas Martin.) 


‘ 
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EXERCISES. 
335 
Prove the formula 
feosz.J4,() 2—2*)dz, 
0 
in which -/, denotes the ordinary Bessel’s function of the zero order. ° 
[Maxime Bécher.] 
336 


Ir Sylvester's difference-product determinant FF, 8) = 
and if we represent this determinant when the 7th column and pth row are 
deleted by J (7, p); evaluate the ratio J (7, p)/J when is infinite. 

[W. Echols.) 
337 

Required the locus of the foot of the perpendicular from the centre of an 
ellipse upon the common chord of the ellipse and circle of curvature. 

[| Artemas Martin.} 
338 

Prove synthetically that the eccentricity of a conic section is equal to the 
sine of the angle which the cutting plane makes with the base of the cone, 
divided by the sine of the angle which an element of the cone makes with the 


base. B. Newson. | 
339 


Snow that the areas of the curves 
y= pert+gqe+r, and y= m2 + pe’ 4+ 
taken between the limits  -- / and  — A, are given by the formule 


Y= + 3 ph*®, and 2 = + p + ma)’, 
respectively. [W. Echols.) 
340 


Fixp the differential equation of plane curves of the third degree free 
from the constants of the general integral equation. {ale Prize Problem.| 


341 


Requirep the form of a beam of uniform strength, supported at its ends, 
the weight of the beam being the only load. [De Volson Wood. | 
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